home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / dlaed6.z / dlaed6
Encoding:
Text File  |  1998-10-30  |  2.8 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAEEEEDDDD6666((((3333FFFF))))                                                          DDDDLLLLAAAAEEEEDDDD6666((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAED6 - compute the positive or negative root (closest to the origin) of
  10.      z(1) z(2) z(3) f(x) = rho + --------- + ---------- + ---------  d(1)-x
  11.      d(2)-x d(3)-x  It is assumed that   if ORGATI = .true
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DLAED6( KNITER, ORGATI, RHO, D, Z, FINIT, TAU, INFO )
  15.  
  16.          LOGICAL        ORGATI
  17.  
  18.          INTEGER        INFO, KNITER
  19.  
  20.          DOUBLE         PRECISION FINIT, RHO, TAU
  21.  
  22.          DOUBLE         PRECISION D( 3 ), Z( 3 )
  23.  
  24. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  25.      DLAED6 computes the positive or negative root (closest to the origin) of
  26.                       z(1)        z(2)        z(3) f(x) =   rho + --------- +
  27.      ---------- + ---------
  28.                      d(1)-x      d(2)-x      d(3)-x
  29.            otherwise it is between d(1) and d(2)
  30.  
  31.      This routine will be called by DLAED4 when necessary. In most cases, the
  32.      root sought is the smallest in magnitude, though it might not be in some
  33.      extremely rare situations.
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      KNITER       (input) INTEGER
  38.                   Refer to DLAED4 for its significance.
  39.  
  40.      ORGATI       (input) LOGICAL
  41.                   If ORGATI is true, the needed root is between d(2) and d(3);
  42.                   otherwise it is between d(1) and d(2).  See DLAED4 for
  43.                   further details.
  44.  
  45.      RHO          (input) DOUBLE PRECISION
  46.                   Refer to the equation f(x) above.
  47.  
  48.      D            (input) DOUBLE PRECISION array, dimension (3)
  49.                   D satisfies d(1) < d(2) < d(3).
  50.  
  51.      Z            (input) DOUBLE PRECISION array, dimension (3)
  52.                   Each of the elements in z must be positive.
  53.  
  54.      FINIT        (input) DOUBLE PRECISION
  55.                   The value of f at 0. It is more accurate than the one
  56.                   evaluated inside this routine (if someone wants to do so).
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAAEEEEDDDD6666((((3333FFFF))))                                                          DDDDLLLLAAAAEEEEDDDD6666((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      TAU          (output) DOUBLE PRECISION
  75.                   The root of the equation f(x).
  76.  
  77.      INFO         (output) INTEGER
  78.                   = 0: successful exit
  79.                   > 0: if INFO = 1, failure to converge
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.